Non— Markovian Quantum Fluctuations and Superradiance Near a Photonic Band 

Edge 
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We discuss a point model for the collective emission of light from N two-level atoms in a photonic 
bandgap material, each with an atomic resonant frequency near the edge of the gap. In the limit of a 
low initial occupation of the excited atomic state, our system is shown to possess novel atomic spectra 
and population statistics. For a high initial excited state population, mean field theory suggests a 
fractionalized inversion and a macroscopic polarization for the atoms in the steady state, both of 
which can be controlled by an external d.c. field. This atomic steady state is accompanied by a 
non-zero expectation value of the electric field operators for field modes located in the vicinity of the 
atoms. The nature of homogeneous broadening near the band edge is shown to differ markedly from 
that in free space due to non-Markovian memory effects in the radiation dynamics. Non-Markovian 
vacuum fluctuations are shown to yield a partially coherent steady state polarization with a random 
phase. In contrast with the steady state of a conventional laser, near a photonic band edge this 
coherence occurs as a consequence of photon localization in the absence of a conventional cavity 
mode. We also introduce a classical stochastic function with the same temporal correlations as the 
electromagnetic reservoir, in order to stochastically simulate the effects of vacuum fluctuations near 
a photonic band edge. 

PACS numbers: 42.50.Fx, 42.50.Lc, 42.70.Qs 



I. INTRODUCTION 

In recent years, photonic bandgap (PBG) structures 
have been shown to lead to the localization of light 
through the carefully engineered interplay between mi- 
croscopic scattering resonances and the coherent inter- 
ference of light from many such scatterers J2L Since 
the initial proposal of photonic bandgaps MM, PBG 
materials exhibiting photon localization have been fab- 
ricated at microwave frequencies Q and more recently, 
large scale two-dimensional PBG systems have been pro- 
duced in the near-infrared The ultimate goal for 
laser applications is a full three-dimensional PBG at op- 
tical frequencies ]t| — |lO|| - A PBG comprises a range of 
frequencies over which linear photon propagation is pro- 
hibited. Therefore, atoms with transition frequencies 
within the gap do not experience the usual fluctuations 
in the electromagnetic vacuum that are responsible for 
spontaneous decay. Instead, a photon-atom bound state 
is formed GJ. Unlike the suppression of spontaneous 
emission from an atom in a high-Q optical microcavity 
fl2|| , the bound photon may tunnel many optical wave- 
lengths away from the atom before being re-absorbed. 
Near a photonic band edge, the photon density of states 
is rapidly-varying, making it dramatically different from 
the w 2 dependence found in free space. This implies that 
the nature of vacuum fluctuations and thus of sponta- 
neous emission near a band edge is radically different 
from that of the exponential decay found in free space 
p3[ . More fundamentally, the correlation time of the 
electromagnetic vacuum fluctuations near a band edge is 
not negligibly small on the time scale of the evolution of 



an atomic system coupled to the electromagnetic field. 
In fact, the reservoir exhibits long-range temporal corre- 
lations, making the temporal distinction between atomic 
system and electromagnetic reservoir unclear. This ren- 
ders the usual Born-Markov approximation scheme in- 
valid for band edge systems. Studies of single atom spon- 
taneous emission near a photonic band edge have 
shown that this non-Markovian system-reservoir inter- 
action gives rise to novel phenomena, such as oscillatory 
behavior and a fractional steady state population for a 
single excited atomic state, as well as vacuum Rabi split- 
ting and a sub-natural linewidth for atomic emission. 

for the collec- 



We consider the Dicke model |16 
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tive emission of light, or superradiance, from N identi- 
cal two-level atoms with a transition frequency near a 
photonic band edge. The study of superradiant emission 
is of interest not only in its own right, but also because 
it provides a valuable paradigm for understanding the 
self-organization and emission properties of a band edge 
laser. Of late, there has been a resurgence of interest in 
superradiance in the context of superradiant lasing ac- 
tion pq ], and due to the experimental realization of a 
true Dicke superradiant system using laser-cooled atoms 
p9| . A low threshold microlaser operating near a pho- 
tonic band edge may exhibit unusual dynamical, spectral 
and statistical properties. We will show that such effects 
are already evident in band edge collective spontaneous 
emission. A preliminary study of band edge superra- 
diance for atoms resonant with the band edge |^0| has 
shown that for an atomic system prepared initially with 
a small collective atomic polarization, a fraction of the su- 
perradiant emission remains in the vicinity of the atoms, 
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and a macroscopic polarization emerges in the collective 
atomic steady state. In addition to this form of sponta- 
neous symmetry breaking, it has been demonstrated that 
superradiant emission can proceed more quickly and with 
greater intensity near a photonic band edge than in free 
space. In the absence of an initial atomic polarization, 
the early stages of superradiance are governed by fluctua- 
tions in the electromagnetic vacuum near the band edge. 
These fluctuations affect the dynamics of collective decay 
and will determine the quantum limit of the linewidth of 
a laser operating near a photonic band edge. 

The organization of the paper is as follows. In Section 
II, we present the quantum Langevin equations for col- 
lective atomic dynamics in band edge superradiance. In 
Section III, we calculate an approximate, analytic solu- 
tion for the equations that describe the N-atom system 
with low initial inversion of the atomic population. We 
show that the atoms can exhibit novel emission spectra 
and a suppression of population fluctuations near a band 
edge. Sections IV and V treat the case of high initial in- 
version. In Section IV, the mean field results of Ref. pc| ] 
are extended to the case of atoms with resonant frequen- 
cies displaced from the band edge. It is shown that the 
phase and amplitude of the collective atomic polarization 
can be controlled by an external field that Stark shifts the 
atomic transition relative to the band edge. The dissi- 
pative effect of dipole dephasing is also included in the 
framework of our non-Markovian system. Section V de- 
scribes superradiant emission under the influence of vac- 
uum fluctuations by exploiting the temporal division of 
superradiance into quantum and semi-classical regimes. 
We find that the system exhibits a macroscopic steady 
state polarization amplitude with a phase precession trig- 
gered by band edge quantum fluctuations. In Section VI, 
we describe a method for generating a classical stochastic 
function that simulates the effect of band edge vacuum 
fluctuations. We show that, for a sufficiently large num- 
ber of atoms, this classical noise ansatz agrees well with 
the more exact simulations of Section V, and may thus be 
useful in the analysis of band edge atom-field dynamics. 
In Appendix A, We give the details of the calculation of 
the electromagnetic reservoir's temporal autocorrelation 
function for different models of the photonic band edge. 
This correlation function is central to determining the 
nature of atomic decay. 



II. EQUATIONS OF MOTION 

We consider a model consisting of N two-level atoms 
with a transition frequency near the band edge coupled 
to the multi-mode radiation field in a PBG material. For 
simplicity, we assume a point interaction, that is, the spa- 
tial extent of the active region of the PBG material is less 
than the wavelength of the emitted radiation. This is of- 
ten referred to as the small sample limit of superradiance 
p7f . We neglect the spatially random resonance dipolc- 



dipole interaction (RDDI) near the band edge, which may 
have a more important impact on atomic dynamics when 
the atomic transition lies deep within the PBG pO| , pl| . 
Nevertheless, our simplified model should provide a good 
qualitative picture of band edge collective emission. For 
an excited atomic state |2) and ground state |1), the in- 
teraction Hamiltonian for our system can be written as 

H = ^TiAxalax + iU^2gx(a x J \ 2 - J 2 iax), (2.1) 



where ax and a x are the radiation field annihilation 
and creation operators respectively; A x = lu x — lo 2 \ 
is the detuning of the radiation mode frequency uj x 
from the atomic transition frequency u) 2 \. g x = 
(u> 2 id 2 i/h)(h/2eoUJ X V) 1 ^ 2 e x -u c i is the atom-field coupling 
constant, where d 2 \\id is the atomic dipole moment vec- 



tor, V is the sample volume, and e> 



ek.< 
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are the two transverse polarization vectors. The Jij 
are collective atomic operators, defined by the relation 

Jij = J2k=i Wfcfc 01 i = !' 2 ' wnere \i) k denotes the 
ith level of the fcth atom. Using the Hamiltonian (|2.l| ), 
we may write the Heisenberg equations of motion for the 
operators of the field modes, a x (t), the atomic inversion, 
J 3 (£) = J22{t) — Ju(t), and the atomic system's collective 
polarization, J\2(t): 



-j;ax{t) = -iA x a x (t) + g x J 12 (t) 
at 



dt 



J 3 (t) = -2^g x J 2l (t)a x {t) + adj. 



dt 



Ji2(t)=J2gxMt)a x (t). 



(2.2) 



(2.3) 



(2.4) 



We may adiabatically eliminate the field operators by 
formally integrating e quat ion ( ^.2j) a nd substituting the 
result into equations ( |2.3| ) and (2.4). The equations of 
motion for the collective atomic operators are then 



■2 / J 21 {t)J 12 {t)G{t-t)dt 
Jo 

■ 2 J21 (*)»?(*) + adj. 



(2.5) 



Mt)= / J 3 {t)J 12 (t)G{t-t)dt + J 3 (t)g(t). (2.6) 



dt 



Here, r](t) = J^x g\ax(0)e 



is a quantum noise opera- 



tor which contains the influence of vacuum fluctuations. 
G(t — t ) is the time delay Green function, or memory 
kernel, describing the electromagnetic reservoir's average 
effect on the time evolution of the system operators. The 
Green function is given by the temporal autocorrelation 
of the reservoir noise operator, 
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G(t-f') = <T,(t)„t (*')) = £ 



9\ e 



-iA x (t-t') 



(2.7) 



We have made use of the fact that (^a\(0)a\(0)j ~ 0, as 

we are dealing with atomic transition frequencies in the 
optical domain [^3|. In essence, G(t — t ) is a measure of 
the reservoir's memory of its previous state on the time 
scale for the evolution of the atomic system. In free space, 
the density of field modes as a function of frequency is 
broad and slowly varying, resulting in a Green function 
that exhibits Markovian behavior, G(t — t) = %8(t — t), 
where 7 is the usual decay rate for spontaneous emis- 
sion p3| . Near a photonic band edge, the density of 
electromagnetic modes varies rapidly with frequency in 
a manner determined by the photon dispersion relation, 
u>k. We show that this results in long range temporal 
correlations in the reservoir which affect the nature of 
the atom-field interaction. 

In order to evaluate G(t — t ) near a band edge, we first 
make the contin uum approximation for the field mode 
sum in equation (2.7): 



G{t-t) 



d 3 k 



2he (2ir) 3 J w k 



, 1 A 2 - 



-i(u k —u>2i)(t—t ) 



(2. 



In this paper, we use an effective mass approximation to 
the full dispersion relation for a photonic crystal. Within 
this approximation, we consider two models for the near- 
band edge dispersion. The details of the calculation of 
G(t — t ) for each model and a discussion of its applicabil- 
ity is given in Appendix A. In an anisotropic dispersion 
model, appropriate to fabricated PBG materials, we as- 
sociate the band edge with a specific point in fc-space, 
k = ko. By preserving the vector character of the dis- 
persion expanded about ko, we account for the fact that, 
as k moves away from ko, both the direction and mag- 
nitude of the band edge wavevector are modified. This 
gives a dispersion relation of the form: 

lu v =w c ±^(k-k ) 2 . (2.9) 

Here, A = 2c 2 /uj gapi where uj gap is the frequency width 
of the gap. The positive (negative) sign indicates that 
LUk is expanded about the upper (lower) edge of the 
PBG, and u> c is the frequency of the corresponding band 
edge. This form of dispersion is valid for a gap width 
ujgap S> c|k — ko|, meaning that the effective mass re- 
lation is most directly applicable to large photonic gaps 
and for wavevectors near the band edge. Furthermore, 
for a large gap and a collection of atoms which are 
nearly resonant with the upper band edge, it is a very 
good approximation to completely neglect the effects 
of the lower photon bands. The band edge density of 
states corresponding to equation (2.E) takes the form 
p(uj) ~ (uj — LUc) 1 / 2 , uj > uj c , characteristic of a three- 
dimensional phase space. The resulting Green function 
for uj c (t — t ) ^> 1 is 



G A (t~t 



^1/2^^/4+5.(4-*')] 
(t-t'f/ 2 



t > t 



(2.10) 



In addition to the anisotropic photon dispersion model, 
it is instructive to consider a simpler isotropic model. 
In this model, we extrapolate the dispersion relation for 
a one-dimensional gap to all three spatial dimensions. 
We thus assume that the Bragg condition is satisfied for 
the same wavevector magnitude for all directions in k- 
space. This yields an effective mass dispersion of the 
form LUk = u c + A(\k\ — |k |) 2 , which associates the band 
edge wavevector with a sphere in fc-space, |k| = k . 
Strictly speaking, an isotropic PBG at finite wavevec- 
tor I ko I does not occur in artificially created, face cen- 
tred cubic photonic crystals. However, a nearly isotropic 
gap near ko = occurs in certain polar crystals with 
polaritonic excitations Effi |. A simple example of such 
a crystal is table salt (NaCl), which has a polariton 
gap in the infrared frequency regime. The band edge 
density of states in the isotropic model has the form 
p(uj) ~ (uj — ujc)^ 1 / 2 , lu > lu c , the square root singularity 
being characteristic of a one dimensional phase space. 
For the Green function we obtain (see Appendix A), 



Gi(t-t) 



^3/2 e _ i [ 7r /4_ <5c ( t _ t ')] 

(t-t'y/ 2 ■ 



t > t 



(2.11) 



In both ( 2. 10] ) and (2.11), S c — OJ21 — ui c is the de- 
tuning of the atomic resonance frequency from the band 
edge, and f3 a is a constant that depends on the dimen- 
sion of the band edge singularity. In particular, for the 



isotropic model, /3 X = cu 21 d\ 1 j\2TitQ'K i l 2 c? while in 
the anisotropic model, (3^ 2 — wfi^ii/^^o^c (An) 3 ^ 2 . 



III. LOW ATOMIC EXCITATION: HARMONIC 
OSCILLATOR MODEL 

In order to understand the effects of band edge vacuum 
fluctuations, we begin by presenting a simplified model 
that permits an analytic solution, and is applicable to 
a system in which only a small fraction of the two-level 
atoms are initially in their excited state. This discus- 
sion demonstrates how light emission near a photonic 
band edge can give rise to novel atomic dynamics, emis- 
sion spectra, and photon number statistics. We write the 
atomic operators in the Schwinger boson representation 
& 



Mt) b\(t)b 2 (t) 

J 3 (t)->#)&2(*) -&!(*)&!(*). 



(3.1) 
(3.2) 



subject to the constraint on the total number of atoms, 
b\(t)bi(t) + b\(t)b 2 (t) = N. The operators b\(t) and b t (t) 
then describe transitions of the system between the ex- 
cited state (i = 2) and the ground state (i = 1). In the 
limit of low atomic excitation, the state |1) has a large 
population at all times, meaning that we can replace the 
inversion operator by the classical value Js(t) « — TV, and 
that bi(t) can be approximated by bx(t) « vJV. In this 
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case, the initially excited two-level atoms behave like a 
simple harmonic oscillator coupled to the non-Markovian 
electromagnetic reservoir. A form of non-Markovian cou- 
pling similar to that of bosons to the electromagnetic field 
occurs in the context of the output coupling of a cold 
atom Bosc condensate from a trapping potential to the 
propagating modes of an atom laser |24|. This mathe- 
matical analogy may lead to deeper insight into both the 
atom laser problem and photonic band edge dynamics. 
In our model, the Heisenberg equations of motion (2.5) 
and (2.6), reduce to 



— b 2 {t) = -N / b 2 (t)G(t-t)dt 
at J 



N V (t). (3.3) 



Using the method of Laplace transforms, we can solve for 
62 (i) and find 

b 2 (t) = B{t)b 2 {Q) - VNj2Mt)^x(0), (3.4) 



where 



B 



(t) = C- 1 {b(s)} , 



B(s) = s + G(s) 



and 



A x (t) = £- 



9\ 



s + iAx 



B(s) 



(3.5) 
(3.6) 

(3.7) 



C^ 1 denotes the inverse Laplace transformation, and 
G(s) is the Laplace transform of the general memory ker- 
nel, G(t — t'). In this section, we consider the case of an 
isotropic b and e dge in the effective mass approximation 
(equation ( 2.1C )), for which G(s) is written as: 



GAs) 



iV/3 3/ V 



■i-rr/4 



\/s — iS c 



(3.8) 



For this isotropic Green functi on, w e denote the inverse 
Laplace transform of equation (3.5) by Bi(t). Bj(t) was 
computed in Ref. |l5| in the context of single atom spon- 
taneous emission, and a detailed mathematical derivation 
may be found therein. Here, it describes the mean or drift 
evolution of our Heisenberg operator b 2 (t). The solution 
has the form 

Bi(t) = 2 cnxxel 31 * 2 ^" 1 + a 2 (x 2 + y 2 )e^ t+lS ' t 
1 - $ 



J2 



(y/Ktfi)] e^ t+ ^\ (3.9) 



where 



(A+ + A_y 



t/4 



x 2 



A. 



T/6-j e -i7r/4 



(3.10) 



(3.11) 



x 3 = (A+e™ /6 - A-e-"/ 6 )^ 3 ^ 4 , 



1 /o\ !/3 



j = 1,2,3. 



(3.12) 



(3.13) 



(3.14) 



$ (x) is the error function, <I> (x) = J* e 4 dt. 

The probability of finding the atoms in the excited 
state is given by (b 2 (t)b 2 (t)^ = \Bi(t)\ 2 , and is plotted 

in Fig. 0. We find that the excited state population ex- 
hibits decay and oscillatory behavior before reaching a 
non-zero steady state value due to photon localization. 
These effects are due to the strong dressing of the atoms 
by the radiation field near a photonic band edge, result- 
ing in dressed atomic states that straddle the band edge. 
Light emission from the dressed state outside the gap re- 
sults in highly non-Markovian decay of the atomic pop- 
ulation, while the dressed state shifted into the gap is 
responsible for the fractional steady state population of 
the excited state. The consequences of this strong atom- 
field interaction are discussed in detail for single atom 
spontaneous emission in Ref. [ p!5[ , and for superradiant 
emission in Sections IV and V of this paper. We note 
that the degree of steady state localization is a sensitive 
function of the detuning, S c , of the atomic resonance from 
the band edge. The decay rate scales as N 2 / 3 (3it for the 
isotropic model. However, there is no evidence for the 
build-up of inter-atomic coherence, as very few of the 
atoms arc in itial ly excited. 

Equation ( |3.3| ) also allows us to calculate the system's 
emission spectrum into the modes lo for an atom with 
resonant frequency uj 2 \ using the relation 

S(UJ) = J e -^-^l)T (r ) ^ (Q^ dj . + C ( , 

~Re{B*[i(u-L02i)}}, (3.15) 



where B(s) is defined in equation (3.E). The spectrum 
for the isotropic model is then 







Nft 



3/2 



LO > UJ C 



(3.16) 



This spectrum is shown in Fig. g, and differs significantly 
from the Lorentzian spectrum for light emission in free 
space. In fact, the emission spectrum is not centered 
about the atomic resonant frequency, which is what one 
would expect for an atom decaying to an unrestricted 
vacuum mode density. We see that for an arbitrary de- 
tuning, 5 C , of lj 2 i from the band edge, the emission spec- 
trum vanishes for frequencies at the band edge and within 
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the gap, u> < lu c . This is consistent with the localization 
of light near the atoms for electromagnetic modes within 
the PBG. As u>2i is detuned further into the gap, spec- 
tral results confirm that a greater fraction of the light is 
localized in the gap dressed state, as the total emission 
intensity out of the decaying dressed state is reduced. 
Conversely, as u>2i is moved out of the gap, the emission 
profile becomes closer to a Lorentzian in form and the 
total emitted intensity increases. The spectral linewidth 
ratio between the isotropic band edge and free space is 
of the order of /Ji/^iV 1 / 3 ), while for an anisotropic band 
edge it is ~ iV/33/7. This corresponds to the fact that col- 
lective emission is much more rapid near an anisotropic 
band edge than in free space, whereas it is slower than 
in free space for the isotropic model. 

It is also instructive to evaluate the quantum fluctua- 
tions in the atomic inversion in the context of the har- 
monic oscillator model. Variances in the atomic popula- 
tion can be written in terms of the Mandel Q-parameter 



Q(t) 



(n*(t)) (n(t)Y 
(n(t)) 



(3.17) 



where n(t) = b\(t)b2(t) is the number operator for the 
occupation of the excited state. Since both the free space 
and PBG solutions in our model can be written in the 
form of equation ( |3.4| ), we can write the Q-parameter in 
the general form 



Q(t) = \B(t)\ 2 Q(0) + Nj2\Mt)f 



(3.18) 



Again, |-B(f)| is the normalized probability of find- 
ing the initially excited fraction of the atoms still in 
the excited state at time t. For an isotropic band 
edge, B{t) = Bi(t) (equation (|3.9[ )), whereas in free 
space, B(t) ~ e~ Njt / 2 , representing the exponential de- 
cay of the excited state population. Using the identity 
NJ2\ \ A x{t)\ 2 = 1 - \ B (t)| , as derived in Appendix B, 
we can write the population fluctuations as 



Q(t) = \B(t)\ 2 [Q(0)-l} + l. 



(3.19) 



For arbitrary initial statistics, atoms in free space de- 
cay to the vacuum state with Q (t) = 1; since the atoms 
decay fully, there are no meaningful atomic statistics in 
the long time limit. Q (t) is plotted in Fig. || for the 
isotropic band edge (S c — 0) for the cases Q(0) = 0, 1, 
and 2. Near the band edge, photon localization prevents 
the atomic system from decaying to the ground state. 
We find instead that the steady state statistics are sensi- 
tive to the statistics of the initial state and to the value 
of S c . A system initially prepared with super-Poissonian 
statistics (Q(0) > 1) experiences a suppression of popula- 
tion fluctuations in the steady state. In a system that is 
initially sub-Poissonian (Q(0) < 1), the fluctuations in- 
crease, but are held below the Poissonian level by photon 



localization. In both cases, the steady state value of the 
atomic population fluctuations is controlled by 8 C . Our 
harmonic oscillator model thus suggests that a PBG sys- 
tem may exhibit novel quantum statistics in the absence 
of a cavity or external fields. It is important to extend 
the analysis of collective emission under the influence of 
vacuum fluctuations to the high excitation (superradi- 
ant) regime. In this case, the two-level nature of the 
atomic operators will become important and will modify 
the quantum statistics from that of the harmonic oscilla- 
tor picture. This generalization is considered in the next 
two sections. 



IV. HIGH ATOMIC EXCITATION: MEAN FIELD 
SOLUTION 

When the atomic system is initially fully or nearly fully 
inverted, we expect inter-atomic coherences, transmitted 
via the atomic polarizations, to have a strong influence 
on emission dynamics. For such high initi al a tomic exci- 
tation, the quantum Langevin equations a nd fl2.6| ), 



paired with the non-Markovian memory kernels ( 2.10| ) or 
( [2.11D , do not possess an obvious analytic solution. More- 
over, conventional perturbation theory applied to these 
equations fails to recapture the influence of the photon- 
atom bound state fiT|| , which plays a crucial role in band 
edge radiation dynamics. However, when the superradi- 
ant system is prepared with a non-zero initial polariza- 
tion (Ji2(0) 0), the average dipole moment dominates 
the incoherent effect of the vacuum fluctuations and the 
subsequent evolution is well-described by a semi-classical 
approximation |17|. In this case, it is possible to factorize 
the atomic operator equations: 



^-(J 3 (t)) = -4Re{(J 21 (t)) f (J 12 (t))G(t-i)dt 



dt 



dt 



(Jl2(t)) 



(4.1) 



(Mt)) y o (Ji 2 (t))G(t-t)dt. (4.2) 



The brackets (O) denote the quantum mechanical aver- 
age of the Heisenberg operator O over the Heisenberg 
picture atom-field state vector, = \vac) ® 1^), where 
\vac) represents the electromagnetic vacuum state, and 
\ip) represents the initial state of the atomic system. 
Clearly in this mean field approach, the quantum noise 
contribution is neglected, as (rj(t)) = 0. Recently, Bay, 
Lambropoulos and Maimer |2q] found that, for a simpler 
Fano profile gap model, the dynamics of superradiant 
emission are affected by the choice of factorization ap- 
plied to the full quantum equations. However, the com- 
plete factorization used here retains the qualitative fea- 
tures and evolution time scales of more elaborate factor- 



ization schemes. Equations (4.1) and (4.2) were solved 
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numerically in reference |20| for an atomic resonance fre- 
quency coincident with the band edge (S c = 0) and a 
small initial collective polarization. The initial collective 
state was assumed to be of the form 

N 

\iP) = \{(V~r\l) + Vl~r\2)) k (4.3) 

fc=l 

with r <C 1, so that initially the atoms are almost fully in- 
verted. In this paper, we extend the previous analysis to 
atomic frequencies detuned from the band edge. Despite 
its neglect of vacuum fluctuations, mean field theory il- 
luminates many of the interesting features of the system. 
The relationship between mean field theory and a more 
complete description including quantum fluctuations is 
discussed in Section V. 

For clarity, we discuss separately the atomic dynamics 
in our isotropic and anisotropic dispersion models. Fig- 
ures |J and ^| show the inversion per atom and the average 
polarization amplitude per atom respectively for various 
values of S c near an isotropic band edge . We see from 
Fig. ^ that a fraction of the superradiant emission re- 
mains localized in the vicinity of the atoms in the steady 
state, due to the Bragg reflection of collective radiative 
emission back to the atoms. This localized light exhibits 
a non-zero expectation value for the field operator, which 
in turn leads to the emergence of a macroscopic polar- 
ization amplitude in the steady state. We further note 
that the decay rate for the upper atomic state is propor- 
tional to iV 2 / 3 . Accordingly, the peak radiation intensity 
is proportional to TV 5 / 3 . This is to be compared with the 
values N and TV 2 for the free space decay rate and peak 
radiation intensity respectively. 

As in single atom spontaneous emission near an 
isotropic band edge [ )15| , the dressing of the atoms by 
their own radiation field causes a splitting of the band 
of collective atomic states such that the collective spec- 
tral density vanishes at the band edge frequency. The 
strongly-dressed atomic states are repelled from the band 
edge, with some levels being pulled into the gap and the 
remaining levels being pushed into the electromagnetic 
continuum outside the PBG. In the long time (steady 
state) limit, the energy contained in the dressed states 
outside the bandgap decays whereas the energy in the 
states inside the gap remains in the vicinity of the emit- 
ting atoms. It is the localized light associated with 
the gap dressed states which sustains the fractionalized 
steady state inversion and non-zero atomic polarization. 
For the isotropic model, this splitting and fractional lo- 
calization persist even when u>2i lies just outside the gap 
(S c > 0) , and the fraction of localized light in the steady 
state increases as u>2i moves towards and enters the gap. 
In the dressed state picture, the self-induced oscillations 
in both the inversion and the polarization which occur 
during radiative emission can be interpreted as being due 
to interference between the dressed states. The oscilla- 
tion frequency is proportional to the frequency splitting 
between the upper and lower collective dressed states. 



This is the analogue of the collective Rabi oscillations of 
N Rydberg atoms in a resonant high-Q cavity (27). From 
Fig. ||, we see that a dressed state outside the band 
gap decays more slowly for atomic resonant frequencies 
deeper inside the gap, causing the collective oscillations 
to persist over longer periods of time. Clearly, this decay 
is non-exponential and highly non-Markovian in nature. 
Fig. confirms that, as required, the polarization ampli- 
tude for large negative values of 5 C is constrained by the 
condition, (J 12 {t)} /N < 1/2. 

In Fig. ^, we plot the phase angle of the collec- 
tive atomic polarization in the isotropic model, 9(t) = 
tan^ 1 {Im (J\2{t)) / Re (Jnif))} ■ Prior to atomic emis- 
sion, this phase angle rotates at a constant rate, and in 
the vicinity of the decay process 9(t) exhibits the effects 
of collective Rabi oscillations. When the emission is com- 
plete, the rate of change of phase angle, 6 if), attains a 
new steady state value, 9(t s ), that depends sensitively on 
the detuning frequency S c . 9{t s ) is a measure of the en- 
ergy difference between the bare atomic state and the lo- 
calized dressed state, Ti {u>2i ~ ^>loc)- Such a polarization 
phase rotation implies that the collective atomic Bloch 
vector of the system exhibits precessional dynamics in 
the steady state. Unlike the conventional precession |2S[] 
of atomic dipoles in an ordinary vacuum driven by an ex- 
ternal laser field, Bloch vector precession in a PBG occurs 
in the absence of an external driving field. Instead, the 
precession is driven by the self-organized state of light 
generated by superradiance, which remains localized near 
the emitting atoms. We see in Fig. ^| that for values of 
5 C such that lo 2 i — oji oc < 0, 9(t s ) is negative, while for 
u>2i —wioc > 0, Q(t s ) is positive, i.e. the phase is rotating 
in the opposite direction. At a detuning corresponding 
to a constant phase in the steady state (9(t s ) — 0), the 
dressed and bare states are of the same energy; this oc- 
curs for a detuning value of S c = — 0.644A r2 / 3 /3i. At this 
value of 5 C , we also find that (J3(t s )) = 0, implying that 
there is no net absorption of light by the atomic system. 
This is, in essence, a collective transparent state [ |28| . 

Collective emission dynamics near an anisotropic band 
edge are pictured in Figs. and |[ For u>2i coincident 
with the band edge or slightly within the gap (S c < 0), 
we again find a fractional atomic inversion in the steady 
state (Fig. Rabi oscillations in the atomic popu- 

lation are much less pronounced than in the isotropic 
model, even for u>2i detuned into the gap. This demon- 
strates that the dressed atomic states outside a physical 
photonic band edge decay much more rapidly than the 
isotropic model would suggest. Furthermore, in contrast 
with the isotropic model, we see that photon localization 
is lost for even a small detuning of u>2i into the continuum 
of field modes outside the band edge. Therefore, while we 
find a macroscopic steady state polarization and preces- 
sional dynamics of the Bloch vector for S c < (Fig. ||), 
for 5 C > the polarization dies away after collective emis- 
sion has taken place. Photon localization from an atomic 
level lying just outside the gap in a three dimensional 
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PBG material may, however, be realized through quan- 
tum interference effects if there is a third atomic level 
lying slightly inside the gap p^) . These results point to 
the greater sensitivity of the atomic dynamics to the more 
realistic anisotropic band edge. Because the isotropic 
model overestimates the momentum space for photons 
satisfying the Bragg condition, photon localization ef- 
fects and vacuum Rabi splitting are exaggerated in the 
isotropic model relative to an artificial photonic crystal. 
In the anisotropic model, the phase space available for 
propagation vanishes as the optical frequency approaches 
the band edge. As a result, vacuum Rabi splitting pushes 
the collective atomic dressed state into a region with a 
larger density of electromagnetic modes. Consequently, 
the decay rate of the atomic inversion is proportional 
to TV 2 near the anisotropic band edge, and the corre- 
sponding peak radiation intensity is proportional to TV 3 . 
Clearly, superradiance near an anisotropic PBG can pro- 
ceed more quickly and can be more intense than in free 
space. As a result, PBG superradiance may enable the 
design of mirrorless, low-threshold microlasers exhibiting 
ultrafast modulation speeds. A more complete analysis 
of lasing near a photonic band edge including pumping 
and dissipative effects will be presented elsewhere. 

From polarization phase and amplitude results, we 
conclude that: (i) Unlike in free space, the atoms near 
a photonic band edge attain a fractionally inverted state 
with constant polarization amplitude and rate of change 
of phase angle. This corresponds to a macroscopic atomic 
coherence in the steady state analogous to that experi- 
enced in a laser. In our case however, "lasing" occurs in 
the band edge continuum rather than into a conventional 
cavity mode, (ii) By varying the value of S c , one may con- 
trol the direction and rate of change of the steady state 
polarization phase angle. This may be realized by apply- 
ing a small external d.c. field to the sample which Stark 
shifts the atomic transition frequency of the atoms. This 
type of control over the collective atomic Bloch vector 
may be of importance in the area of information storage 
and optical memory devices |30|j3l[ |. 

The above analysis makes it clear that collective spon- 
taneous emission dynamics in a PBG are significantly 
different from those in free space. In a real PBG ma- 
terial, the dephasing of atomic dipoles due to inter- 
atomic collisions or phonon-atom interactions may also 
have a significant effect on the evolution of our system 
over a large range of temperatures. In the free space 
Markov approach, dipole dephasing is described by a phc- 
nomenological polarization decay constant fl32| . Since the 
Markov approximation does not apply near a band edge, 
one cannot account for dephasing by si mply adding a phc- 
nomenological decay term to equation (4.2). However, we 
expect that the atomic resonant frequency will experience 
random Stark shifts due to atom-atom or atom-phonon 
interactions. This effect can be included in the descrip- 
tion of our system by adding a variation A to the de- 
tuning frequency 5 C at e ach time s tep in a computational 
simulation of equations (4.1) and (4.2). A is chosen to be 



a Gaussian random number with zero mean. The width 
of the Gaussian distribution is determined by the mag- 
nitude of the random Stark effect. Such a simulation in 
free space would include a random A only in the equation 
for the atomic polarization. This is because the slowly 
varying photon density of states seen by the atoms at 
the frequency u>2i + A does not change significantly with 
typical homogeneous line broadening effects. In contrast, 
we have seen that near a photonic band edge, slight vari- 
ations in 5 C may drastically change the atomic inversion. 
Therefore we include A in both system equations. In 
Fig. |], we plot the evolution of the collective inversion 
and polarization under the simulated collision broaden- 
ing described above. The random Stark shifts lead to the 
loss of macroscopic polarization and the loss of atomic in- 
version in the long time limit. The latter effect can be 
understood by noting that the random frequency shifts 
are symmetrically distributed about the mean resonant 
frequency. Frequency shifts into the gap promote photon 
localization, while those away from the gap cause further 
decay of the atomic inversion. Over time, the net result 
is that the frequency shifts away from the gap encourage 
the decay of the atomic population. This is true even in 
atomic systems for which the mean resonant frequency 
lies within the gap. From the above considerations, it is 
clear that dephasing is a significant perturbation on pho- 
ton localization near a photonic band edge. In particular, 
dephasing determines the threshold external pumping re- 
quired to achieve atomic inversion in a band edge laser. 
It also facilitates the emission of laser light from the pho- 
tonic crystal. 

Although a superradiant system can be prepared in 
a coh eren t initial state of the type described by equa- 
tion (4.3) 1 28 1, collective emission is typically initiated 
by spontaneous emission, a random, incoherent process. 
Over time, spontaneous emission leads to the build-up 
of macroscopic coherence in the sample. The effect of 
vacuum fluctuations is then of considerable importance 
in the full description of superradiance, both from a fun- 
damental point of view, and for potential device appli- 
cations, such as the recently proposed superradiant laser 
]i~8"| . In the next section, we present a more detailed de- 
scription of PBG superradiance that takes into account 
the role of quantum fluctuations. 



V. BAND EDGE SUPERRADIANCE AND 
QUANTUM FLUCTUATIONS 

In order to describe the evolution of the superradiant 
system's collective Bloch vector under the influence of 
quantum fluctuations, we consider atomic operator cor- 
relation functions of the form |33l 



9 Pq = {(Jl2f(J2l) q ) 



(5.1) 



Here the operators are evaluated at equal times. As in 
free space, we expect vacuum fluctuations to drive the 
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system from its unstable initial state with all atoms in- 
verted to a new stable equilibrium state. Such fluctu- 
ations are particularly relevant prior to the build-up of 
macroscopic atomic polarization. Indeed, they provide 
the trigger for superradiant emission. In the early-time, 
i nve rted r egim e, we may set 7 3 (t) — J 3 (0) in equations 
O) and (ph, giving " 



Ju(t) = / dfJ 3 (0)Jx 2 (t )G(t -t) + J 3 (0)v(t). 



dt 



(5.2) 



The resulting equation remains non-linear, and involves 
products of atomic and reservoir operators. We may 
simplify expressions containing operators in this inverted 
regime by considering operator averages over only the 
atomic Hilbert space. For an arbitrary Heisenberg oper- 
ator 0(t), we denote the atomic expectation value for an 
initial fully inverted state |7) by (0) A = (I\0\I). We 
denote by the set {| A)} a complete set of 2 N normalized 
basis vectors for the atomic Hilbert space including |7), 
such that (A|7) = S\j, where 5 a ,/3 is the Kronecker delta 
function. Clearly, (7 | J 3 (0) | A) = N5j,\. Since J 3 (0) 
acts as a source term for J\2{t) in equation (5.2), we also 
have the property (7 | J12 (t) | A) = for A 7^ I in the 
inverted regime. This can be shown by considering the 
equation of motion for (7| J\2{t) |A): 

dt"£(I\M0)\fi) ( M | J 12 (f) |A) G(t-t') 
(7| J 3 (0)|A)ry(t) 



= N dt'(I\J u (t')\\)G(t-t')+N5 ItX ri(t), (5.3) 
Jo 

where fj, labels a complete set of atomic states. This 
integro-differential equation satisfies the initial condition 
(7 1 Ji2(0) |A) = 0, since Ji2(0) acts as a raising operator 
on the fully inverted bra vector (7|. For A / /, the 
source term in ( |5.3| ) is also absent, leading to the so- 
lution (7| Jvz{t) I A) = 0. Using this property, we may 
replace the atomic average over products of atomic op- 
erators with products of atomic averages, provided that 
7 3 (t) — 73(0). For example, 

(J12J21) = ( vac \ ® i 1 1 Jl2 I M) (Ml J21 1 7) <g> \vac) 



((Jl2) A (J2l) A ) ] 



(5.4) 



Here, (O) R = (vac\ O \vac) denotes an expectation value 
over the reservoir variables. For an arbitrary moment 
g pq , we have 



^ = ((7i 2 )^(7 21 )«). 



(5.5) 



We note that such a factorization is valid only for an 
antinormal ordering of the polarization operators, since 
(7| J2i(i) I A) does not vanish in general. 



Taking the atomic expectation value of equation (|j) , 
we obtain 

j t (Ji2(t)) A = N J* (Mt')) A G(t - t)dt + Nr,{t). 

(5.6) 

This is a linear equation that has lost its operator charac- 
ter over the atomic variables but not over the electromag- 
netic reservoir, as evidenced by the pr esen ce of the quan- 
tum noise operator, rj(t). Equation ( |5.6| ) can be solved 
by the method of Laplace transforms. The solution has 
the form, 

(Mt)) A - D(t) (J 12 (0)) A + Nj2Cx(t)a x (0), (5.7) 

A 

where, 

D(t) = C- 1 [b{s)) , (5.8) 



D(s) = s - NG(s) 



and 



C x (t) = £- 



<-l\ 



s + iAx 



D(s) 



(5.9) 



(5.10) 



Again, C~ l denotes the inverse Laplace transform. The 
Laplace transformation of the memory kernel fo r a n 
isotropic band edge, G/(s), is given in equation (|3.g|). 
Despite the fact t hat (^2(0))^ = 0, we retain the first 
term in equation (5.7) for later notational convenience 

The early-time quantum fluctuations in a superradiant 
system prevent us from predicting a priori the evolution 
of any single experimental realization of the atoms. In- 
stead, we can only determine the probability of a partic- 
ular trajectory of the collective atomic Bloch vector. In 
order to obtain the statistics of a band edge superradi- 
ant pulse, we first determine the statistics of the collec- 
tive Bloch vector for a set of identically-prepared systems 
after each has passed through the early time regime gov- 
erned by vacuum fluctuations. The relevant time scale 
will be referred to as the quantum to semi-classical evolu- 
tion crossover time, t — tg. Our approach is to calculate 
the phase and amplitude distributions of the polariza- 
tion at the crossover time quantum mechanically. The 
subsequent (t > to) evolution of the ensemble is the n ob- 
tained by solving the semi-classical equations J4.l| ) and 
(4.2) using the polarization distribution function at to. 
In other words, the distribution of values of (Ji2(to)) ob- 
tained from the early time quantum fluctuations provide 
the initial conditions for subsequent, semi-classical evo- 
lution. In order to implement this approach, we must 
first identify £0 for our system jl7],[34]]. One expects such 
a transition to occur in the high atomic inversion regime, 
(J3(t)) ~ N. It is natural to define to such that for t > to 
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the expectation value of the commutator of the system 
operators J 2 i(t) and Ji 2 (t) becomes very small compared 
to the expectation value of their product jl7| . This gives 
the condition, 

{Jn{t)Jxi{t)) » ([J21W, Jl2(t)]) , i > to. (5.11) 

Evaluating the above commutator, we have 
([J 2 i(t), Jvi(t)]) = (Mt )), w hich is equal to TV for full 
atomic inversion. From (5.S) and (5.7), we find that 



(J2i{t)J 12 (t)) = (J 3 (t) + Ji 2 (t) J 2 i(t)> 



iV 



l + i\r£|C A (t)| 5 



= A|£(t)| 2 , 



The last equality is obtained by use of the identity 
N Y,x \Cx(t)\ 2 = \D(t)\ 2 - 1, as derived in Appendix B. 
In free space, |D(t)| 2 = e Nlt , giving the crossover time, 
ff ree ~ 1/7V7. One can solve for the crossover time near 
a band edge, t BBG , computationally. In the isotropic 
model, for S c = we find that t PBG ~ 1.24/iV 2 / 3 /3 1 . The 
crossover time maintains this l/iV 2 / 3 /3i dependence for 
u>2i displaced from the band edge. The corresponding 
time scale for the anisotropic gap is 1/N 2 (3 3 . The build- 
up of a macroscopic polarization then occurs more slowly 
near an isotropic and more quickly near an anisotropic 
band edge than in free space. 

Using a semi-classical approach, we may write the 
value of the polarization at any time t > to in terms 
of an amplitude k and a phase <j>, (J\2(t)) Cl = J(k, <j>, t). 
The superscript CI refers to the fact that the expectation 
value ( ) Cl is taken in the semi-classical regime t > to. 
We define P{n)dn as the probability of finding the ampli- 
tude between k and K+dn, and Q(<j))d4> as the probability 
of finding the phase between <fi and </> + d<f>. We may then 
write the moments of the macroscopic polarization dis- 
tribution as 



(Mt)7 (J2i(t)) 9 ) 



CI 



dK d(j)P{K)Q{(j))[J{K,(j) 1 t)] P J*(K,<f),t) 



(5.13) 

For t = to, we assume that the polarization has the form 
J(k, (f>, t) — ne 1 ^ , giving for the moments 

((JMY (J2i(t )) q f l 

dK { d<j>P{K)Q{<j))K p+q e i ^- q) ' t > . (5.14) 



The quantum analogue, ( ) , of (5.14) can be written 
in the f orm of equation (5J3) eval uate d at t = to. Substi- 
tuting (5.7) and its adjoint into (5.5) yields: 



((Jl2(to)) P (J2l(to)) 9 ) 



q\Q 



JSfP+i 



E C A(io)al(0) 



' R 

(5.15) 



As the reservoir expectation value is taken over the the 
operators a\ which satisfy a Gaussian probability distri- 
bution, Wick's theorem [[t3| is applied in order to reduce 
the operator averages of products of field operators to 
averages over products of pairs of field operators. We 
then have 

<(Ji 2 (to)) P (J 2 i(to)) 9 ) Q * S„N*p\ \D(t Q )\ 2p . (5.16) 



(5.12) This expression has corrections of order N p 1 , mean- 



ing that it is asym ptotically valid for large N. Equat- 
ing ( 5.14 ) and ( 5.16 ), we solve for the distributions P(k) 
and Q(4>) to obtain the desired initial polarization dis- 
tribution for the semi-classical superradiance equations. 
The early time distributions for free space and the band 
edge differ only in the form of the function D(t), as the 
above analysis makes no other distinction between the 
two cases. Thus in the band edge system, as in free space, 
the entire effect of the early time atomic evolution can 
be recaptured using the distribution of initial conditions 
given at t = to. The phase of the polarization is given by 
the relation 



2tt 



(5.17) 



This shows that Q(<f) is uniformly distributed between 
and 2n. The initial polarization amplitude distribution 
is found from the relation 



dnP (k) k 2p 



N\D(t )\< 



(5.18) 



The result is a Gaussian distribution of width N\D(to)\' 
centered at zero, 



P(«) = 



1 



yJ*N\D{to)f 



: exp 



N\D{to)\' 



(5.19) 



It has been shown via density matrix methods [ p"7| that 
in free space one may choose the crossover time any- 
where in the inverted regime, the simplest choice being 
to = 0. This is due to the absence of temporal corre- 
lations of the reservoir for t ^ t . Figure |l0| shows the 
ensemble-averaged collective emission in free space and 
at an isotropic band edge (S c — 0) for N=100 atoms. 
Both the free space and band edge systems are shown for 
two choices of initial polarization distribution. The solid 
lines corresp ond t o the choice of to = in the amplitude 
distribution ( [5.1£ ) for both free space and the band edge. 



The dashed lines correspond to the choice to = to ree and 







to = t^ BG for the free space and b and edge systems re- 
spectively. As per equation ( [3. 17 ), the initial phase of 
the polarization in all cases is chosen from a uniform 
random distribution. In free space, the two choices of 
initial conditions yield indistinguishable atomic dynam- 
ics. This verifies that the choice of to is unimportant in 
free space, so long as it is chosen in the inverted regime. 
Near a photonic band edge, we see that the choice of 
to affects the later evolution of the system. In particu- 
lar, it affects the onset time for collective emission. It 
is clear from these simulations that the details of the 
non-Markovian evolution in the quantum regime play a 
crucial role in the subsequent semi-classical evolution of 
the band edge superradiance. The long-range tempo- 
ral correlations of the reservoir require that we treat the 
vacuum fluctuations explicitly throughout the quantum 
evolution of the system. A similar picture holds in the 
case of an anisotropic PBG material. In our anisotropic 
model, memory of the initial state is expressed through 
the Green function ( |2.11 ). In this case, superradiance is 
also highly sensitive to early stage quantum fluctuations. 

Since ensemble averages of atomic observables are ex- 
perimentally measurable quantities, we consider these in 
some detail. We use the notation ( ) ens to denote an 
ensemble-averaged quantum expectation value. For il- 
lustration, we focus on the S c — and zero dephasing 
case for a system of 100 atoms in the isotropic effec- 
tive mass model. The extension to non-zero detuning 
and finite dipole dephasing follows from the discussion 
of Section IV. From Fig. [y], it is evident that the en- 
semble exhibits a fractional population inversion in the 
steady state. The steady state value of (J3(t)) ens for a 
given atomic detuning is unchanged from the mean field 
result, (J3(t s )). Since the steady state is determined by 
the atom-field coupling strength, and not by the dynam- 
ics of the system, it is insensitive to initial conditions. 
Fluctuations in the excited state atomic population may 
be expressed in terms of the delay time for the onset 
of superradiant emission, defined as the time at which 
the system is exactly half-excited, i.e. (J3) = 0. Vac- 
uum fluctuations result in a distribution of delay times 
for the ensemble, asymmetrically centered about a peak 
value, as pictured in Fig. O. The delay time distribu- 
tion is qualitatively similar to that obtained in free space 
p4j . However, the width of the distribution scales with 
the relevant time scale for the isotropic and anisotropic 
gaps, showing that, near a photonic band edge, atomic 
population fluctuations during light emission can be re- 
duced from their free space value. Because of the varia- 
tion in initial conditions, the Rabi oscillations in (Js(t)) 
for the isotropic gap are much less pronounced than in 
mean field simulations. The differences in emission times 
due to fluctuations cause the ensemble average inversion 
to smear out these oscillations. Therefore, one can no 
longer directly relate the amplitude and period of the os- 
cillations to the energies of the collective dressed states. 

More striking is the nature of the ensemble's collective 
polarization under the influence of vacuum fluctuations. 



Figures |l3| a - d show the evolution of the polarization 
distrib ution from the initial distribution given by equa- 
tions ( |5.17 ) and ( 5.19 ) to the steady state distribution. 
Initially, the distribution is sharply peaked about zero. In 
the decay region, the polarization amplitude is broadly 
distributed and has a random phase. This behavior is 
reminiscent of the fluctuations of the order parameter in 
the vicinity of a phase transition. In the steady state, the 
polarization amplitude collapses to a very well-defined 
non-zero value. This amplitude is again accompanied by 
a random phase that is uniformly distributed between 
and 2tt. We may interpret our steady state result in 
the following manner: A fraction of the photons emit- 
ted near the photonic band edge remain localized in the 
vicinity of the atoms, causing both the atomic dipoles and 
the electromagnetic field to self-organize into a cooper- 
ative steady state. However, vacuum fluctuations cause 
this cooperative quantum state to have a random phase, 
resulting in a zero ensemble average polarization ampli- 
tude, |(Ji2(t)) e „ s | = 0, as shown in Fig. |ll|. Measure- 
ments of the degree of first and second order coherence 
of the electromagnetic field in a band edge superradiance 
experiment would provide a probe of the nature of this 
self-organized state of photons and atoms near a band 
edge. Such measurements would however have to be done 
by indirect means, as the localized radiation field in a 
PBG does not propagate out of the photonic crystal. We 
further note that this state - well-defined in amplitude 
but with random phase - is similar to the steady state of 
a conventional laser |35| with a well-defined electric field 
and random phase diffusion. 



VI. SIMULATED QUANTUM NOISE NEAR A 
BAND EDGE 

We have shown that the statistical properties of a band 
edge superradiant system can be determined because the 
collective behavior of the constituent atoms leads to a 
semi-classical system evolution, triggered by early-time 
quantum fluctuations. However, a seamless quantum de- 
scription of band edge quantum optical systems is ex- 
tremely difficult to obtain, due to the non-Markovian 
nature of the atom-field interaction. As a first step, we 
introduce a method by which to simulate their evolu- 
tion computationally and include the effects of quantum 
fluctuations. Unlike the semi-classical simulations of sec- 
tion IV which neglected the effect of the quantum noise 
operator, as (r](t)) — (rf(t)) — 0, we propose to replace 
(r](t)} in our semi-classical equations by a complex classi- 
cal stochastic function with the same mean and two-time 
correlation function as its quantum counterpart. This 
noise function then simulates the quantum noise in our 
system throughout the entire system evolution. We may 
test the validity of our simulated noise ansatz for band 
edge superradiance by comparing the results obtained to 
those calculated in Section V. 
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The classical noise function required to simulate quan- 
tum noise near a photonic band edge involves a real 
stochastic function £ (t) possessing the underlying tempo- 
ral autocorrelation of our non-Markovian quantum noise 
operator, rj{t). In the effect ive m ass a pprox imation, this 
means that (see equations ( 2.1C ) and ( 2.11 )), 



(6.4) 



mat 



i 



(6.1) 



where again a = 1 and 3 for isotropic and anisotropic 
band edges respectively. Problems in band edge atom- 
field dynamics, such as the present superradiant prob- 
lem, often involve non-linear equations under the influ- 
ence of colored quantum noise. It is interesting to note 
that nonlinear problems involving classical colored noise 
are of considerable interest in classical statistical physics. 
In particular, methods fo r co mputationally generating 
noise satisfying equation (6.1) have been developed in 
the context of problems in surface growth and poly- 
mer ordering |36[ . In what follows, we use the method 
first introduced by Rice |57j and elaborated on by Bil- 
lah and Shinozuka |38| for the generation of such colored 
noise. For noise with a power spectrum P(u), defined 
as the Fourier transform of the autocorrelation function 
(£(*)£(*'))> their algorithm gives, 

N 

£(t) ^ 2^[PK)A w ] 1/2 cosK< + $„), 

71=1 

n = l,2,...,JV, (6.2) 

with equality obtained for N — ► oo. Here, u> n = nAu, 
Aw = w max /7V, and u) max is a cutoff frequency above 
which the power spectrum can be neglected. Each is 
a random phase uniformly distributed in the range [0, 27r] 
. By use of a particular set of random phases to 
generate the noise values at each time step, we obtain a 
single " experimental" realization of the quantum noise in 
our system. Since we cannot predict a priori the specific 
form of the quantum fluctuations in a particular exper- 
iment, we again average over many realizations of the 
superradiant system, each governed by a different 
in order to obtain distributions and ensemble averages of 
relevant quantities. We note that equation (6.2) clearly 
gives {£(t)) = 0, as desired, since the random $„ cause 



the ensemble to average to zero. To show that (3.2) also 
gives the correct autocorrelation function, we write: 

( mm ) ens 

k I 

xcos(w,t' + $i) ) ens 

= (2A^^[PK)/ 3 h)f /2 

k I 

x {cos (uo k t - (Jit' + 3> fc - $z) 

+ cos ((j k t + uit' + $ fc + $i)}) ens (6.3) 



In ( |3.3j ), only the k — I components, in which the random 
phases and $; cancel each oth er, survive the ensemble 
average. All other terms in ( |6.3|) vanish in the ensemble 
average. As N — > oo, d6.4| ) becomes the Fourier trans- 
form of P((j), which equals (£(*)£(*')) ens ■ Studies have 
shown that for values of N as small as 1000, the desired 
autocorrelation may be obtained with as little as 5% error 
p8| , making this a computationally feasible technique. 
Furthermore, unlike other methods for the generation of 
stochastic functions (see Refs. f|(|), the present method 
computes the desired function, using only a uni- 

form random distribution of phases <&k as input, rather 
than requiring the computation of a Gaussian stochas- 
tic function as an intermediate step. This decreases the 
likelihood of spurious correlations between our random 
numbers. Figure |l4| shows the two-time correlation of 
£(t) for a = 1, for an ensemble of 2000 realizations of 
the noise function generated by the algorithm of equa- 
tion (6.2). In this calculation and in the simulations de- 
scribed below, we chose a power spectrum P(u>) — yjji 
in order to mimic the colored vacuum near an isotropic 
band edg e. W e see good agreement with the correlation 
function (6J). The agreement between our simulations 
and the exact correlation function can be significantly 
improved by enlarging the size of the ensemble, at the ex- 
pense of increased computation time for atom-field sim- 
ulations. 

The ensemble {£(£)} is used t o sim ulat e th e effect of 
vacuum fluctuations in equations (2.5) and (2.6). Written 
in terms of the dimensionless time variable r = N 2 / 3 (3\t, 
these equations for the isotropic band edge become 



^( J 3(T)) 



(Jai(r)> r( J2l(T,) )^(r-rV 



-4i?e{e-" /4 - _ , 

V 71 " JO VT - T 

{J2i{T))e- l ^/ % - s ^ 



<(r) } 



(6.5) 



dr V 71 " J s/t - t 

(6.6) 



(J 3 (r)) e-W- 6 '^ 



C(r), 



with similar equations for the anisotropic gap. For both 
models, the noise term scales as 1/viV; this is the same 
dependence of the noise term o n p article number exhib- 
ited in free space |}3|. In Fig. |lj, we show the average 
inversion for an ensemble containing 2000 realizations of 
£(t) for N = 1000 and N = 10000 atoms. We find 
that our stochastic simulation scheme gives physical re- 
sults only for systems of N > 500 atoms. The stochastic 
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simulations show good agreement with the atomic inver- 
sion obtained by the method of Section V. Other system 
properties, such as the ensemble-averaged polarization 
and the delay time distribution calculated by the present 
method also agree well with the quantum calculations of 
the previous section. This suggests that our stochastic 
approach may be a valuable tool in the analysis of band 
edge atom-field dynamics. 

VII. CONCLUSIONS 

In this paper, we have treated the collective sponta- 
neous emission of two-level atoms near a photonic band 
edge. An analytic calculation of the atomic operator dy- 
namics in the case of low atomic excitation was given. 
The results demonstrate novel atomic emission spectra 
and show the possibility of reducing atomic population 
fluctuations. This in turn suggests that fluctuations in 
photon number are likewise suppressed for light localized 
near the atoms. This raises the interesting question of 
whether squeezed light |3j|, antibunched photons fl40| , 
and other forms of non-classical light may be generated 
in a simple manner from band edge atom-field systems. 
For an initially inverted system prepared with a small 
macroscopic polarization, a mean field factorization was 
applied to the atomic quantum Langevin equations, giv- 
ing a semi-classical system evolution. We found that 
the atoms exhibit fractional population trapping and a 
macroscopic polarization in the steady state. Collective 
Rabi oscillations of the atomic population were found, 
and were attributed to the interference of strongly- 
dressed atom-photon states that are repelled from the 
band edge, both into and out of the gap. The degree of 
photon localization, the polarization amplitude, and the 
phase angle of the polarization in the steady state are 
all sensitive functions of the detuning of the atomic res- 
onant frequency from the band edge. The steady state 
atomic properties can thus be controlled by applying a 
d.c. Stark shift to the atomic resonant frequency. 

The effect of quantum fluctuations for high initial ex- 
citation of the atoms was included by distinguishing 
regimes of quantum and semi-classical collective atomic 
evolution. We found that the early time quantum evolu- 
tion must be treated in detail, due to the non-Markovian 
electromagnetic reservoir correlations near a band edge. 
This is in contrast with free space, where the atomic sys- 
tem's evolution is insensitive to the treatment of the full 
temporal evolution of the early, quantum regime. Frac- 
tional localization of light was shown to persist under 
the influence of vacuum fluctuations. The atomic polar- 
ization exhibits a non-zero amplitude with a randomly 
distributed phase in the steady state. This is much like 
the steady state of a conventional laser. Here, such las- 
ing characteristics are due only to the Bragg scatter- 
ing of photons back to the atoms; there is neither ex- 
ternal pumping nor a laser cavity in our system. The 



time scales for all dynamical processes, such as collective 
emission and the buildup of collective atomic polariza- 
tion are strongly modified from their free space values 
due to the singular photon density of states near a pho- 
tonic band edge. For an isotropic band edge, the time 
scales as iV 2 / 3 /?i, while in the more realistic anisotropic 
model, time scales as N 2 j3^. As a result, collective emis- 
sion phenomena can occur more rapidly near a band edge 
than in free space. Throughout our calculations, we have 
employed an effective mass approximation to the band 
edge dispersion. For materials with a very small PBG, 
it may be important to include the effects of both band 
edges. These issues are raised in Appendix A. 

We have demonstrated that band edge superradiance 
possesses many of the self-organization and coherence 
properties of a conventional laser. Furthermore, we have 
shown the possibility for the generation of novel emis- 
sion spectra and photon statistics. These results sug- 
gest that a laser operating near a photonic band edge 
may possess unusual spectral and stastical properties, as 
well as a low input power lasing threshold due to the 
fractional inversion of the atoms in the steady state. It 
may further be possible to produce a PBG laser in a 
bulk material without recourse to a defect-induced cav- 
ity mode. Lending credence to this hypothesis, recent 
observations EJ and theoretical studies E3] of lasing 
from a multiply-scattering random medium with gain 
have demonstrated that one may obtain light with the 
properties of a laser field in the absence of a cavity. We 
are currently investigating the properties of a band edge 
laser, including the effects of pumping and dissipation, 
such as the dipolar dephasing modeled in Section IV. A 
full description of the statistics of a band edge laser field 
will likely require a non-perturbative master equation (or 
its equivalent) that exhibits the non-Markovian nature 
of the electromagnetic reservoir. Techniques for treat- 
ing the atom-field interaction in the absence of the Born 
and Markov approximations have recently been obtained 
|^3|,^|. However, these methods cannot account for the 
van Hove singularity in the density of states encountered 
at a photonic band edge, nor are they directly applicable 
to externally driven atomic systems. In the absence of 
a more exact approach, the band edge stochastic noise 
function described in Section VI may be used to recover 
many of the results of a full quantum treatment. 

Finally, we note that the steady state atom-field prop- 
erties described here are a result of the effect of radiation 
localized in the vicinity of the active two-level atoms. 
This leads to the question of how to pump energy into 
these states, which lie within the forbidden photonic gap. 
One possibility is to couple energy into and/or out of the 
system through a third atomic level whose transition en- 
ergy lies outside the gap . There is also the possibility 
of transmitting light into the gap through high intensity 
ultrashort pulses that locally distort the nonlinear dielec- 
tric constant of the material and thus allow the propa- 
gation of light in the form of solitary waves within the 
forbidden frequency range [Q. Such issues must be ad- 
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dressed in order to fully exploit the very rich possibilities 
of quantum optical processes near a photonic band edge. 
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APPENDIX A: CALCULATION OF THE 
MEMORY KERNEL 

We first present the calculation of the memory kernel 
for the isotropic model in the effective mass approxima- 
tion, Gi(t — t'). Starting from equation ( pTg| ) and the 
isotropic dispersion relation near the upper band edge, 
ojk = oj c + A(\\c\ — |k |) 2 , G[{t — t') can be expressed as 



4He TT 



A e - iA (k-k ) 2 (t-t ) 



dfl I dk 



A(k-k ) 2 
(Al) 



Again, S c — &21 — ui c is the detuning of the atomic res- 
onant frequency from the band edge. A = mc/h is a 
cutoff in the photon wavevector above the electron rest 
mass. Photons of energy higher than the electron rest 
mass probe the relativistic structure of the electron wave 
packets of our resonant atoms |^(| . Because the isotropic 
model associates the band edge with a sphere in fc-space, 
there is no angular dependence in the expansion of lo^ 
about the band edge. We may thus se par ate out the an- 
gular integration over solid angle fi in (Al). We may also 
make a stationary phase approximation to the integral, 
as the non-exponential part of the integrand will only 
contribute significantly to the integral for k ~ fco- The 
resulting integral is 



G 



f ( t _ t ') ~ *O e «e(t-f ) / dk e -iA{k-k nt-t')_ 

4fte 7r uj c Jk 

(A2) 



As A is a large number, we extend the range of integration 
to infinity in order to obtain a simple analytic expression 
forGj(t-tf) @: 

Cr/(t — t) = 



12^e 7r 3 / 2 c 3 Vt~~F 

3/2 e -*[V*-«-(*-0] 



(A3) 



Because the relevant frequencies in equation ( |A3| ) 
are roughly of the same order of magnitude near a 
band edge, we may re-write the prefactor as (3^ 2 ~ 

w 2i 2 ^2i/12fi.eo7r 3 / 2 c 3 , in agreement with the value given 
in Section II. We emphasize that the stationary phase 
method yields the correct asymptotic behavior for the 
memory kernel for large \t — t'\. At short times, the in- 
tegral must be evaluated more precisely using the full 
photon dispersion relation, as discussed below. 

For an anisotropic band gap model, we must ac- 
count for the variation in the magnitude of the band 
edge wavevector as k is rotated throughout the Bril- 
louin zone. We associate the gap with a specific point 
in fc-space that satisfies the Bragg condition, k = k . 
In the effective mass approximation, the dispersion rela- 
tion is expanded to second order in k about this point, 
= oj c ±A{k.~ ko) 2 . Making the substitution q = k— ko 
and performing the angular integration, G a {t — t ) is ex- 
pressed as 



G A (t-t = — -e 



21 "21 iSJt-i) I j„9 e 



A 2 -iAq 2 (t-t ) 

dq— t^T- ( A4 ) 

u c + Aq 2 



Extending the wavevector integration to infinity, the 
Green function is |47[| 



G A (t - 1 ) 



, , 2 rl 2 
^21»21 

8fte 7r 2 



,iS c (t-t ) 



H i Wc (t-t') 

2V A 



iio c (t — t') 



}■ 

(A5) 

$ (x) is the probability integral, $ (x) = e~ l dt. 

For u) c [t — t ) 3> 1, a condition satisfied for all but the 
t — > t limit (as uj c ~ 10 15 s _1 for optical transitions), 
taking the asymptotic expansion of <fr(x) to second order 
gives 



G A (t-t 



, i 2 rl 2 
w 21«21 



j[ir/4+S c (t-t')\ 



8he {irA) 3 / 2 Lj c (t-t') 3 / 2 ' 

oj c (t - t) > 1. 

(A6) 



As t 

G A (t-t) = 



t — > 0+, ( A5) reduces to 

/ i 2 rl 2 



8he ir 2 A 3 / 2 



-TTy/uTc 

VjJ c \t — t J 



t-t 



-0+. 
(A7) 



G A (t — t ) possesses a weak (square root) singularity at 
t = t . This is an integrable singularity and can thus be 
treated numerically p8| . 

The effective mass dispersion relation used in the eval- 
uation of G(t — t ) is, strictly speaking, valid only near 
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the photonic band edge, as it fails to give the required 
linear photon dispersion relation for large \k — ko\ (far 
away from the gap). Therefore, the integration of th e ef- 
fective mass dispersion for large wavevector in (|Al[) and 
( A4) introduces a spurious contribution to G(t — t ). This 
difficulty may be overcome for an isotropic gap model by 
using a dispersion relation that has the correct wavevec- 
tor dependence for all k. The simplest model dispersion 
with the correct form is 

LOk/c = ^Jkl + 7 2 + sgn{k - fc )\/(fc - fc ) 2 + 7 2 - (A8) 

The double-valued nature of ajfc at k is made explicit 
by the function sgn(k — ko) = 1 for k > ko, and 
sgn(k — ko) = —1 for k < kg. This gives a gap of 
width lu 



with G(t — t') defined as in equation (2.7). We may re- 
write this double integral in the form: 



dt' / dt"B(t")B*(t')G(t' -t") 
o Jo 

+ / dt' [ dt"B{t")B*{t')G{t' -t") 
Jo Jf 

h+h, (B3) 



where I± and li are the first and second double integrals 
respectively. By changing the order of the integrations 
in I2, we obtain 

t" 



gap 



27c, centered about the midgap frequency I 2 = / dt" dt 1 B(t")B* (t')G(t' - t") = If. (B4) 

Jo Jo 



^0 = c y/ko + 7 2 - Note that equation ( |A8| ) gives the cor- 
rect linear dependence in k for both large positive and 
negative k, and gives the effective ma ss d ispersion for 
k ~ fco- Like the effective mass model, (A8) gives a sin- 
gular density of states at the band edges, oj c = luq ± cy. 
The full dispersion relation allows us to evaluate the in- 
fluence of both band edges for arbitrary gap width and 
atomic resonant frequency. Preliminary numerical cal- 
culations show a stronger reservoir memory effect than 
demonstrated in the effective mass model for the isotropic 
band edge. This may have a significant effect on theoreti- 
cal predictions regarding the atom-field interaction in the 
vicinity of a PBG. A further simplification has been made 
in the anisotropic model, in that we have not included the 
dependence of Wk on the symmetry of a specific photonic 
crystal. In a real three dimensional PBG material, the 
Bragg condition is satisfied for different values of k as the 
wavevector changes direction in fc-space. This directional 
dependence may lead to a much stronger dependence of 
the localization of light on the detuning of 1U21 away from 
the band edge. The impact on the atom-field interaction 
in a PBG of both the full isotropic dispersion model and 
more realistic dispersions for three dimensional photonic 
crystals will be treated elsewhere. 



Therefore, ^ A | j 4a(0T — {/1}, and we need only ex- 
plicitly eva luat e I\. The Laplace transform of B(t), B(s) 
(equation (^6)), is equivalent to the Laplace transform 
of the equation 



dt 



B{t) =-N dt'G{t - t')B(t'). 



(B5) 



Substituting ( p35| ) into 1\ and its complex conjugate, we 
obtain 



2/i 



1 / dt i± m >)f = ± 
NJo dt' 1 y n 



N 



\B(0)\ 2 -\B(t)\ 2 



(B6) 



as l\ is real. Substituting the initial condition |i?(0)| z 



1 into (B6) gives the result quoted in Section III. 



APPENDIX B: EVALUATION OF J2x 2 

We outline the evaluation of ^ A |A,\(£)| 2 , used to ob- 
tain the low e xcitat ion population fluctuations in Section 
III, equation (|3.19|). A similar procedure is used to ar- 
rive at equation ( 5.12 ) in Section V. Starting from the 
Laplace transform A\(s) (equation ( |3.7| )), we may use 
the properties of a convolution of Laplace transforms in 
order to write 



' J3f4-'\„ — 'Ait 



Ax(t)=g x / dt'B{t')e 
Jo 



(Bl) 



Therefore, we have 



J2\A x (t)\ 2 = f dt' f dt"B{t")B*(t')G(t' -t"), (B2) 
Jo Jo 
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FIG. 1. Normalized population of the excited atomic state 
near an isotropic photonic band edge for low initial atomic 
excitation. Various values of the detuning, 5 C = u>2i — w c , 
of the atomic resonant frequency u>2i from a band edge at 
frequency lu c are shown. Dashed line, 5 C = —.5; solid line, 
S c = 0; dotted line, S c = .5. S c is measured in units of N 2 ^ 3 (3i. 

FIG. 2. Collective atomic emission spectrum <S(w) (arbi- 
trary units) near an isotropic band edge for low initial atomic 
excitation. Various values of the detuning, 5 C = wi\ — ijJ c , of 
the atomic resonant frequency o»2i from an isotropic photonic 
band edge at frequency lu c are shown. Dotted line, S c — — 1; 
dashed line, S c 5 = 0; solid line, 5 C = 1. S c is measured in 
units of iV 2/3 /3i. 

FIG. 3. Fluctuations in the excited state atomic popula- 
tion as measured by the Mandel parameter, 
Q(t) = ((n 2 (i)) - (n(t)) 2 )/(n(t)), for low initial excitation 
for an atomic resonant frequency tuned to an isotropic pho- 
tonic band edge, 8 C = 0. Dashed line, Q(0) = 2; solid line, 
Q(0) = 0. Long-short dashed line denotes fluctuations for 
Poissonian population variance, Q(0) = 1. 

FIG. 4. Mean field solution for the atomic inversion, 
(Js(t)) /N, near an isotropic photonic band edge, starting 
with an infinitesimal initial polarization, r = 10~ 5 . Various 
values of the detuning, S c = LO21 — co c , of the atomic resonant 
frequency LO21 from a band edge at frequency lj c are shown, 
(a) 8 C = 1; (b) 8 C = .5; (c) 8 C = 0; (d) 8 C = -.5; (e) 8 C = -1. 
S c is measured in units of N 2 ' A f3\. 
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FIG. 5. Mean field solution for the atomic polarization am- 
plitude, | («/i2 (*)) | /N, near an isotropic photonic band edge, 
starting with an infinitesimal initial polarization, r = ICP 5 . 
Various values of the detuning, 8 C = u>2i — <^ c , of the atomic 
resonant frequency W21 from a band edge at frequency u) c are 
shown, (a) S c = 1; (b) S c = .5; (c) 5 C = 0; (d) 8 C = — .5; (e) 
S c = — 1. S c is measured in units of N 2 ^ 3 (3\. 

FIG. 6. Mean field solution for the phase angle of the 
atomic polarization, 6(t), near an isotropic photonic band 
edge, starting with an infinitesimal initial polarization, 
r — ICP 5 . Various values of the detuning, 8 C = W21 — oj c , 
of the atomic resonant frequency W21 from a band edge at fre- 
quency uj c are shown, (a) 8 C — .5; (b) 8 C — 0; (c) S c — —.75; 
(d) 8c = — 1. 8 C is measured in units of N 2 ^ 3 (i\. 

FIG. 7. Mean field solution for the atomic inversion, 
(J3(t)} /N, near an anisotropic photonic band edge, starting 
with an infinitesimal initial polarization, r = ICP 6 . Various 
values of the detuning, 8 C = W21 — <^ c , of the atomic resonant 
frequency u)2i from a band edge at frequency u) c are shown. 
Dashed line, 8 C = .1; solid line, 8 C — 0; dotted line, 8 C = —.3. 
S c is measured in units of N 2 (3'j,. 

FIG. 8. Mean field solution for the atomic polarization am- 
plitude, I (J12 (i)} ! /N, near an anisotropic photonic band edge, 
starting with an infinitesimal initial polarization, r — ICP 6 . 
Various values of the detuning, 8 C = W21 — <^c, of the atomic 
resonant frequency L021 from a band edge at frequency lo c are 
shown. Dashed line, 8 C = .1; solid line, 8 C = 0; dotted line, 
8 C — — .3. 8 C is measured in units of N 2 ^. 



FIG. 12. Distribution of delay times for a system of 100 
atoms at an isotropic band edge (8 C = 0) for 2000 realizations 
of the superradiant system. 

FIG. 13. Atomic polarization distribution for a system of 
100 atoms at an isotropic band edge (8 C = 0), subject to 
quantum fluctuations at early times. 5000 realizations of the 
superradiant system, (a) t = to PBG ; (b) t = 5; (c) t = 11; (d) 
steady state, t in units of 1/N 2 '' 3 (3\. 

FIG. 14. Solid line: ensemble averaged autocorrelation 
function, (£,{t)£,(t')) ena , of the classical colored noise function 
£(t) corresponding to vacuum fluctuations near an isotropic 
band edge. The dashed line is a plot of the exact au- 
tocorrelation function in the effective mass approximation, 
(t-tT 1/2 . 

FIG. 15. Comparison of the ensemble averaged atomic in- 
version, (J r 3(i)) ens /N, at an isotropic band edge (8 C = 0) as 
calculated by the methods of Sections V and VI. 2000 realiza- 
tions of the superradiant system. Dashed line and long-short 
dashed line: inversion calculated by the method of Section V 
for N — 1000 and 10000 atoms respectively. Solid line and 
dotted line: inversion calculated using the stochastic function 
of Section VI for N = 1000 and 10000 atoms respectively. 



FIG. 9. Mean field solution for the atomic inversion (solid 
line) and polarization amplitude (dashed line) under the influ- 
ence of collision broadening for an atomic resonant frequency 
at an isotropic photonic band edge, S c = 0. The system is 
given an infinitesimal initial polarization, r = 10 _J . The sim- 
ulated stark shift is a Gaussian random distribution with zero 
mean and standard deviation .5N 2 ^ 3 /3i. 



FIG. 10. Atomic inversion for superradiance driven by vac- 
uum fluctuations in free space and for an atomic resonant fre- 
quency tuned to an isotropic photonic band edge (8 C = 0). 
Solid lines: result for initial polarization distribution at t — 
for each system; dashed lines: result for initial polarization 
distribution at t = to for each system. 



FIG. 11. Ensemble-averaged atomic 

inversion, {Js(t)} ens /N, and atomic polarization amplitude, 
\(JMt)) en3 \ /N (dot-dashed line), for a system of N = 100 
atoms near an isotropic photonic band edge. The ensemble 
average is taken over 2000 initial polarization values. Inver- 
sion: long dashed curve, 8 C = —.5; solid line, 8 C = 0; short 
dashed line, 8 C = .5. S c in units of N 2 ^ 3 P\. 
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